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Abstract 

In our work on the development of model-independent data analysis methods for de- 
termining ratios between different couplings/cross sections of Weakly Interacting Massive 
Particles (WIMPs) by using measured recoil energies from direct Dark Matter detection 
experiments directly, it was assumed that the analyzed data sets are background-free, i.e., 
all events are WIMP signals. In this article, as a more realistic study, we take into account 
a fraction of possible residue background events, which pass all discrimination criteria and 
then mix with other real WIMP-induced events in our data sets. 

Our simulations show that, assuming that the spin-dependent (SD) WIMP-nucleus 
interaction dominates over the spin-independent (SI) one, the maximal acceptable fraction 
of residue background events in the analyzed data sets for determining the ratio of the SD 
WIMP coupling on neutrons to that on protons is ~ 20% - 40%; whereas considering a gen- 
eral combination of the SI and SD WIMP interactions, the maximal acceptable background 
ratio for determining the ratio between two SD WIMP couplings as well as the ratios of 
the SD cross section on protons (neutrons) to the SI one is ~ 10% - 20%. Moreover, by 
considering different forms of background spectrum, we find that only background events 
in the lowest energy ranges could affect the reconstructions (significantly); those in high 
energy ranges would almost not change the reconstructed ratios or only very slightly. 



1 Introduction 



Currently, direct Dark Matter detection experiments searching for Weakly Interacting Massive 
Particles (WIMPs) are one of the promising methods for understanding the nature of Dark 
Matter (DM) and identifying them among new particles produced at colliders as well as recon- 
structing the (sub)structure of our Galactic halo [1, 2, 3, 4]. To this aim, model-independent 
methods for determining the WIMP mass [5, 6], the spin-independent (SI) WIMP coupling on 
nucleons [7, 8] as well as ratios between different WIMP couplings/cross sections [9, 10] from 
direct detection experiments have been developed. 

These methods built basically on the work on the reconstruction of the (moments of the) one- 
dimensional velocity distribution function of halo WIMPs by using experimental data (measured 
recoil energies) directly [11]. The spectrum of recoil energy is proportional to an integral over the 
one-dimensional WIMP velocity distribution, fi{v), where v is the absolute value of the WIMP 
velocity in the laboratory frame. Since this integral is in fact just the minus-first moment of the 
velocity distribution function, which can be estimated from experimental data directly [11, 6], 
by assuming that the spin-dependent (SD) WIMP-nucleus interaction dominates over the spin- 
independent one, an expression for determining the ratio between the SD WIMP coupling on 
neutrons and that on protons has been derived [9, 10]. Meanwhile, for a general combination 
of the SI and SD WIMP-nucleus interactions, by using detector materials with and without 
spin sensitivities on protons and/or on neutrons, a second expression for determining the ratio 
between two SD WIMP-nucleon couplings as well as two expressions for determining ratios of 
the SD WIMP-proton(neutron) cross section to the SI one have also been derived [9, 10]. It was 
found that, by combining experimental data sets with different target nuclei, the ratios between 
different WIMP couplings/cross sections can be determined without making any assumption 
about the velocity distribution of halo WIMPs nor prior knowledge about their mass [9, 10]. 

In the work on the development of these model-independent data analysis procedures for 
extracting information on WIMP couplings/cross sections from direct detection experiments, it 
was assumed that the analyzed data sets arc background-free, i.e., all events are WIMP signals. 
Active background discrimination techniques should make this condition possible. For example, 
the ratio of the ionization to recoil energy, the so-called "ionization yield" , used in the CDMS- 
II experiment provides an event-by-event rejection of electron recoil events to be better than 
10~^ misidentification [12]. By combining the "phonon pulse timing parameter", the rejection 
ability of the misidentified electron recoils (most of them are "surface events" with sufficiently 
reduced ionization energies) can be improved to be < 10~® [12]. Moreover, as demonstrated 
by the CRESST collaboration [13], by means of inserting a scintillating foil, which causes some 
additional scintillation hght for events induced by a-decay of ^^°Po and thus shifts the pulse 
shapes of these events faster than pulses induced by WIMP interactions in the crystal, the pulse 
shape discrimination (PSD) technique can then easily distinguish WIMP-induced nuclear recoils 
from those induced by backgrounds^. 

However, as the most important issue in all underground experiments, the signal identification 
ability and possible residue background events which pass all discrimination criteria and then mix 
with other real WIMP-induced events in analyzed data sets should also be considered. Therefore, 
in this article, as a more realistic study, we follow our works on the effects of residue background 
events in direct Dark Matter detection experiments [17, 18, 19] and want to study how well we 
could determine the ratios of WIMP-nucleon couphngs/cross sections model-independently by 
using "impure" data sets and how "dirty" these data sets could be to be still useful. 

^For more details about background discrimination techniques and status in currently running and projected 
direct detection experiments see e.g., Refs. [14, 15, 16]. 



The remainder of this article is organized as follows. In Sec. 2 I review briefly the model- 
independent methods for determining ratios between different WIMP-nucleon couplings/cross 
sections by using experimental data sets directly. In Sec. 3 the effects of residue background 
events in the analyzed data sets on the measured energy spectrum will be discussed. In Sees. 4 
and 5 I show numerical results of the reconstruction of ratios of WIMP-nucleon couplings/cross 
sections by using mixed data sets with different fractions of residue background events based 
on Monte Carlo simulations. I conclude in Sec. 6. Some technical details will be given in an 
appendix. 



2 Methods for determining ratios of WIMP-nucleon cou- 
plings/cross sections 

In this section 1 review briefly the model-independent methods for determining the ratio of the 
SD WIMP coupling on neutrons to that on protons as well as the ratio between the SD and 
SI WIMP-proton cross sections^. For more detailed illustrations and discussions about these 
precedures see [10]. 



2.1 Event rate for elastic WIMP— nucleus scattering 



Considering the SI and SD WIMP-nucleus interactions together, the basic expression for the 
differential event rate for elastic WIMP-nucleus scattering can be given as [3, 10]: 
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Here R is the direct detection event rate, i.e., the number of events per unit time and unit mass 
of detector material, Q is the energy deposited in the detector, po is the WIMP density near 
the Earth, ctq^^'^^^ are the SI/SD total cross sections ignoring the form factor suppression and 
-^(si,SD)(Q) are the elastic nuclear form factors for the SI/SD WIMP interactions, respectively, 
fi{v) is the one-dimensional velocity distribution function of the WIMPs impinging on the 
detector, v is the absolute value of the WIMP velocity in the laboratory frame. The reduced 
mass ttIt n is defined by 
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where is the WIMP mass and m-^ that of the target nucleus. Finally, Vmin is the minimal 
incoming velocity of incident WIMPs that can deposit the energy Q in the detector: 



with the transformation constant 



a 



(3) 



(4) 



and Vmax is the maximal WIMP velocity in the Earth's reference frame, which is related to the 
escape velocity from our Galaxy at the position of the Solar system, Vesc ^ 600 km/s. 



^In this section I consider only the case with protons, but aU formulae given in Section 2.3 can be modified 
straightforwardly to the case with neutrons. 



Isotope 


Z 


J 


(5p) 


(5n) 


-(.5p)/(5n) 


(5n)/(5p) 


Natural abundance (%) 




9 


1/2 


0.441 


-0.109 


4.05 


-0.25 


100 


23Na 


11 


3/2 


0.248 


0.020 


-12.40 


0.08 


100 




17 


3/2 


-0.059 


-0.011 


-5.36 


0.19 


76 




17 


3/2 


-0.058 


0.050 


1.16 


-0.86 


24 


73 Gc 


32 


9/2 


0.030 


0.378 


-0.08 


12.6 


7.8 / 86 (HDMS) [23] 


127j 


53 


5/2 


0.309 


0.075 


-4.12 


0.24 


100 


i29xe 


54 


1/2 


0.028 


0.359 


-0.08 


12.8 


26 


i3iXe 


54 


3/2 


-0.009 


-0.227 


-0.04 


25.2 


21 



Table 1: List of the relevant spin values of the nuclei used for simulations presented in this 
paper. More details can be found in e.g., Refs. [3, 20, 21, 22]. 

Through e.g., squark and Higgs exchanges with quarks, WIMPs could have a "scalar" inter- 
action with nuclei^. The SI scalar WIMP-nucleus cross section can be expressed as [3, 4] 



Here 
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is the SI WIMP cross section on protons, /p(n) are the effective xxpp(nn) four-point couplings, 
A is the atomic mass number of the target nucleus, and mr,p is the reduced mass of the WIMP 
mass m-^ and the proton mass m^. Note that I have used here the theoretical prediction for the 
lightest supersymmetric neutralino (and for all WIMPs which interact primarily through Higgs 
exchange) that the scalar couplings are approximately the same on protons and on neutrons: 
/n ~ /p; the tiny mass difference between a proton and a neutron has also been neglected. 

On the other hand, through e.g., squark and Z boson exchanges with quarks, WIMPs could 
also couple to the spin of target nuclei, an "axial-vector" (spin-spin) interaction. The SD 
WIMP-nucleus cross section can be expressed as [3, 4]: 



(5'p)ap + (5'n)a„ . (7) 



Here is the Fermi constant, J is the total spin of the target nucleus, (5'(p^n)) are the ex- 
pectation values of the proton and neutron group spins^, and a(p,n) are the effective SD WIMP 



•^Besides of the scalar interaction, WIMPs could also have a "vector" interaction with nuclei [3, 4]. However, 
for Majorana WIMPs (x = x)) e.g., the lightest neutralino in supersymmetric models, there is no such vector 
interaction. 

"^Notc that detailed nuclear spin structure calculations show that not only unpaired nucleons contribute to 
the total cross section, the even group of nucleons has sometimes also a non-negligible spin (see Table 1 and e.g., 
data given in Refs. [3, 20, 21]). 



couplings on protons and on neutrons. Since for a proton or a neutron J — ^ and (Sp) or 
(Sn) — |, the SD WIMP cross section on protons or on neutrons can be given as 

<^xfp,n) = (— ) mj(p_„) |a(p,n) 1^ . (8) 

As shown in Eq. (5), due to the coherence effect with the entire nucleus, the cross section for 
SI scalar WIMP-nucleus interaction scales approximately as the square of the atomic mass of 
the target nucleus. Hence, in most supersymmetric models, the SI cross section for nuclei with 
A > 30 dominates over the SD one [3, 4]. However, as discussed in Refs. [24, 25, 26], in Universal 
Extra Dimension (UED) models, the SD WIMP interaction with nucleus is less suppressed and 
could be compatible or even larger than the SI one. 



2.2 Only a dominant SD WIMP— nucleus cross section 

Consider at first the case that the SD WIMP-nucleus interaction strongly dominates over the SI 
one and thus neglect the first SI term, (Jq^F|j((5), in the bracket on the right-hand side of Eq. (1). 
By using a time-averaged recoil spectrum, and assuming that no directional information exists, 
the normalized one-dimensional velocity distribution function of halo WlMPs, fi{v), has been 
solved analytically [11] and, consequently, its generalized moments can be estimated by [11, 6]^ 
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Here v{Q) = a^/Q^ Q(min,max) are the experimental minimal and maximal cut-off energies of the 
data set, respectively, 

KQmin) ^ (^) (10) 

\ Cxpt, (3 = Qmin 

is an estimated value of the measured recoil spectrum {dR/dQ)cxpt {before normalized by an 
experimental exposure, S) ai Q = Qmin, and In{Qmm, Qmax) can be estimated through the sum: 

A^tot Q{n-l)/2 

-^n(Qmin) Qmax) Tp2( r) \ ' {^^) 

where the sum runs over all events in the data set that satisfy Qa € [<5min, <5max] and A^tot is 
the number of such events. Then, since the integral on the right-hand side of Eq. (1) is just 
the minus-first generalized moment of the velocity distribution function, (f^^), which can be 
estimated by Eq. (9), by setting Q = Qmm and using the definition (4) of a, one can obtain 
straightforwardly that 
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Now, in order to eliminate po here, we combine two experimental data sets with different target 
nuclei, X and Y. By substituting the expression (7) for (7q° into Eq. (12) and using again the 



^Here we have implicitly assumed that Qmax is so large that a term 2(5max^^^^''(<9max)/^'^(Qmax) is negligible. 



definition (4) of a for both target nuclei, the ratio between two SD WIMP-nucleon couplings 
has been solved analytically as [9, 10]^ 
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Here I have used the following relation [6]: 

T^n,Y 



nj^O. 



(13) 



ax 
ay 



and defined 

T^J,n,X = 

with^ 



J 



X 



Jx + llTZnX 



1/2 



(14) 



(15) 



1 



Fx{Qmm,x) 



o,x 



(16) 



and 



T^n.X — 



r,^{n+l)/2 
^Vmin.X 



rx{Qmm,x) / F^{Qmm,x) + (^ + 1)/„,X 



2<3mki,X^x(<3min,x)/-^l(<5min,x) + Io,X 



l/n 



(17) 



TZj,n,Y-, TZa,Y: and Tln,Y Can be defined analogously^; -F(x,y)(Q) are the form factors of the 
nucleus X and Y, T(x,y)(<5min,(x,y)) refer to the counting rates for the target X and Y at 
the respective lowest recoil energies included in the analysis, and S(^x,y) ai'e the experimental 
exposures with the target X and Y. Note that, firstly, Eq. (13) can be used once positive signals 
are observed in two (or more) experiments; information on the local WIMP density po, on the 
velocity distribution function of incident WIMPs, fi{v), as well as on the WIMP mass are 
not necessary. Secondly, because the couplings in Eq. (7) are squared, we have two solutions for 
On/op here; if exact "theory" values for 7^j,n,(x,y) are taken, these solutions coincide for 
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for TZj,n,x = , 
for 7^J,„,y = , 



(18) 



which depends only on properties of two used target nuclei (see Table 1). Moreover, it can be 
found from Eq. (13) that one of these two solutions has a pole at the middle of two intersections, 
which depends simply on the signs of {S^)x and {S^)y- since 7lj^n,x and 7?.j„y are always 



^Note that, although the constraints on two SD WIMP-nucleon couplings have conventionally been shown in 
the flp — On plane, I will always use the On/ctp ratio in this article. 

^Notc that TZcy.(x.Y) ^^nd 7?.„ (x.y) defined here as well as the estimator for /„ given in Eq. (11) can be used 
for either the SI or the SD case with a corresponding form factor. However, since we consider here only the SD 
interaction, F^{Q) needed for using Eqs. (12), (16) and (17) should be substituted by form factors for the SD 
cross section. 

^Hereafter, without special remark all notations defined for the target X can be defined analogously for the 
target Y and occasionally for the target Z. 



positive, if both {Sn)x and (5'n)y are positive or negative, the "— (minus)" solution (on/op)^'^ 
will diverge and the "+ (plus)" solution (an/ap)^^„ will be the "inner" solution; in contrast, if 
the signs of {Sn)x and (S'n)r are opposite, the "— (minus)" solution will be the "inner" solution. 

On the other hand, it has been found [10] that, in order to reduce the statistical uncertainty 
on (an/ap)|°„ estimated by Eq. (13),^. one can practically use the estimate of the counting rate, 
instead of at the experimental minimal cut-off energy, at the shifted point Qs,i (from the central 
point of the first bin, Qi) defined by 



Qs,i = Qi + ^ In 



sinh(A;i6i/2) 



kibi/2 



(19) 



where ki is the logarithmic slope of the reconstructed recoil spectrum in the first Q— bin and bi 
is the bin width. Then, according to Eq. (A9), the measured recoil spectrum at Q = Qs,i can be 
estimated by 



(20) 



where A^^i is the event number in the first bin. 

As shown in Ref. [10], the statistical uncertainties on (an/flp)±°n estimated with different n 
(namely with different moments of the WIMP velocity distribution function) with r(x,y) {Qs,i,(x,y)) 
are clearly reduced and, interestingly, almost equal. Therefore, since 
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one would need practically only events in the lowest energy ranges for estimating a^/a.^. And, 
consequently, one has to estimate the values of form factors only aX. Q = Qs,i, and the zero 
momentum transfer approximation F^((5 ^ 0)) ~ 1 can then be used. 

2.3 Combination of the SI and SD cross sections 



Now I consider the case with a non-negligible SI WIMP-nucleus cross section. At first, by 
combining Eqs. (5), (7), and (8), we can find 
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Then the expression (1) for the differential event rate can be rewritten as 
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^It is true with non-negligible experimental threshold energies [10]. Later we will see that, with negligible 
threshold energies, (on/ap)±^„ estimated with r(^x,Y){Qs,i,{x,Y)) could be a little bit larger. 



Set Q — Qmin- One can find straightforwardly that, for this general case, Eq. (12) becomes to 
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where /„ should be estimated by Eq. (11) with the replacement of F'^{Q) by F''^{Q) = Fgi{Q) + 

By combining two targets X and Y and using the relation (14) for ax /ay with n = —1, the 
ratio of the SD WlMP-proton cross section to the SI one has been solved analytically as [10] 
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where m(x,y) oc A(^x,y) has been assumed, Cp,(x,y) have been defined in Eq. (23), and 
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As the estimator (13) for a^/a-p, one can use Eq. (26) to estimate cr^p/'^xp ^'^^hout a prior 
knowledge of the WIMP mass m-^. Moreover, since Cp^(x,y) depend only on the nature of the 
detector materials, cr^p/'^xp practically only a function of 7^m,(x,y), i-e., the counting rate at 
the experimental minimal cut-off energies, which can be estimated by using events in the lowest 
available energy ranges. 

Meanwhile, for the general combination of the SI and SD WIMP-nucleus cross sections, the 
On/op ratio appearing in Eq. (23) has been solved analytically by introducing a third nucleus 
with only an SI sensitivity: {Sp)z — {Sn)z — 0, i.e., Cp^z = as [10] 
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Here I have defined 
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Note that, firstly, (an/op)^^ and Cp,(x,y) given in Eqs. (28), (29a), and (29b) are functions 
of only r(x,y,z)(Qmin,{x,y,z)) (or ri^x,Y,z){Qs,i,{x,y,z))), which can be estimated with events in the 



lowest energy ranges. Secondly, while the decision of the inner solution of (an/ap)±^ri depends 
on the signs of (S'n)x and (5'n)y, the decision with (an/flp)±^^° depends not only on the signs 
of -Sn/p.x = {Sn)x/ (Sp) X and Sn/p,y = ('S'n)y/ ('S'p)y, but also on the order of the two targets. For 
e.g., a Ge + CI combination, since Sn/p/SQe = 12.6 > Sn/p,37ci = —0.86, one should use the upper 
expression in the second line of Eq. (28), and since Sn/p Tsge and Sn/p,37ci have the opposite signs, 
the "— (minus)" solution of this expression (or the "+ (plus)" solution of the expression in the 
first line) is the inner solution. In contrast, for the F + I combination used in our simulations, 
since s^/p = —0.247 < Sn/p.i^^i = 0.243 and since s^/^wp and s^/p i27j have the opposite signs, 
the "— (minus)" solution of the lower expression in the second line of Eq. (28) (or the "— 
(minus)" solution of the expression in the first line) is then the inner solution. 

Furthermore, in order to reduce the statistical uncertainty, one can choose at first a nucleus 
with only an SI sensitivity as the second target: (<S'p)y = (-S'n)y = 0, i.e., Cp_y = 0. The 
expression in Eq. (26) can thus be reduced to 
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Then we choose a nucleus with (much) larger proton (or neutron) group spin as the first target: 
{Sp)x ^ ('S'n)x — 0, in order to eliminate the On/op dependence of Cp^x given in Eq. (23):^° 
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3 Effects of residue background events 

In this section I first show some numerical results of the energy spectrum of WIMP recoil signals 
mixed with a few background events. For generating WIMP-induced signals, we use the shifted 
Maxwellian velocity distribution [2, 3, 11]: 
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with vo ^ 220 km/s and Ve — 1.05 vq, which are the Sun's orbital velocity and the Earth's 
velocity in the Galactic frame^^, respectively; the maximal cut-off of the velocity distribution 
function has been set as Vmax = 700 km/s. The commonly used elastic nuclear form factor for 
the SI WIMP-nucleus cross section [27, 3, 4]: 
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as well as the thin-shell form factor for the SD WIMP cross section [2, 28, 10]: 

j^iqRi) , for qRi < 2.55 or qR^ > 4.5 , 

const. ~ 0.047 , for 2.55 < qRi < 4.5 



(36) 



Analogously, we can define 
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and choose {Sn)x ^ ('S'p)x — to eliminate its a^/a^ dependence. 

^^The time dependence of the Earth's velocity will be ignored in our simulations. 



will also be used^^. Meanwhile, in order to check the need of a prior knowledge about an (exact) 
form of the residue background spectrum, two forms for the background spectrum have been 
considered. The simplest choice is a constant spectrum: 



More reahstically, we use the target-dependent exponential form introduced in Ref. [17] for the 
residue background spectrum: 



Here Q is the recoil energy, A is the atomic mass number of the target nucleus. The power index 
of A, 0.6, is an empirical constant, which has been chosen so that the exponential background 
spectrum is somehow similar to, but still different from the expected recoil spectrum of the 
target nucleus; otherwise, there is in practice no difference between the WIMP scattering and 
background spectra. Note that, among different possible choices, we use in our simulations 
the atomic mass number A as the simplest, unique characteristic parameter in the general 
analytic form (38) for defining the residue background spectrum for different target nuclei. 
However, it does not mean that the (superposition of the real) background spectra would depend 
simply /primarily on A or on the mass of the target nucleus, m^. In other words, it is practically 
equivalent to use expression (38) or {dR/dQ)\jg^cx = e"'^^^^'^ directly for a ^^Ge target. 

Note also that, firstly, as argued in Ref. [17], two forms of background spectrum given above 
are rather naive; however, since we consider here only a few residue background events induced 
by perhaps two or more different sources, which pass all discrimination criteria, and then mix 
with other WIMP-induced events in our data sets of 0(50) total events, exact forms of different 
background spectra are actually not very important and these two spectra, in particular, the 
exponential one, should practically not be unrealistic^^. Secondly, as demonstrated in Ref. [10] 
and reviewed in the previous section, the model-independent data analysis procedures for deter- 
mining ratios between different WIMP-nucleon couplings/cross sections require only measured 
recoil energies (induced mostly by WIMPs and occasionally by background sources) from di- 
rect detection experiments. Therefore, for applying these methods to future real data, a prior 
knowledge about (different) background source(s) is not required at all. 

Moreover, for our numerical simulations presented here as well as in the next two sections, 
the actual numbers of signal and background events in each simulated experiment are Poisson- 
distributed around their expectation values independently; and the total event number recorded 
in one experiment is then the sum of these two numbers. Additionally, we assumed that all 
experimental systematic uncertainties as well as the uncertainty on the measurement of the 
recoil energy could be ignored. The energy resolution of most existing detectors is so good that 
its error can be neglected compared to the statistical uncertainty for the foreseeable future with 
pretty few events. 

In Figs. 1 I show measured energy spectra (solid red histograms) for a ^^Ge target with six 
different WIMP masses: 10, 25, 50, 100, 250, and 500 GeV based on Monte Carlo simulations. 
The dotted blue curves are the elastic WIMP-nucleus scattering spectra, whereas the dashed 

^^Other commonly used analytic forms for the one-dimensional WIMP velocity distribution as well as for the 
elastic nuclear form factor for the SI WIMP-nucleus cross section can be found in Refs. [11, 29]. 

^^Other (more realistic) forms for background spectrum (perhaps also for some specified targets/experiments) 
can be tested on the AMIDAS website [30, 31]. 
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Figure 1: Measured energy spectra (solid red histograms) for a ^^Ge target with six different 
WIMP masses: 10, 25, 50, 100, 250, and 500 GeV. The dotted blue curves are the elastic WIMP- 
nucleus scattering spectra, whereas the dashed green curves are the exponential background 
spectra normalized to fit to the chosen background ratio, which has been set as 20% here. The 
experimental threshold energies have been assumed to be negligible and the maximal cut-off 
energies are set as 100 keV. The background windows have been assumed to be the same as the 
experimental possible energy ranges. 5,000 experiments with 500 total events on average in each 
experiment have been simulated. See the text for further details (plots from Ref. [17]). 
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Figure 2: As in Figs. 1, except that a ^^F target has been used, the input WIMP mass has 
been fixed as 100 GeV, and four different background ratios: no background (top left), 10% (top 
right), 20% (bottom left), and 40% (bottom right) are shown here. 

green curves are the exponential background spectra given in Eq. (38), which have been normal- 
ized so that the ratios of the areas under these background spectra to those under the (dotted 
blue) WIMP scattering spectra are equal to the background-signal ratio in the whole data sets 
(e.g., 20% backgrounds to 80% signals shown in Figs. 1). The experimental threshold energies 
have been assumed to be negligible and the maximal cut-off energies are set as 100 keV. The 
background windows (the possible energy ranges in which residue background events exist) have 
been assumed to be the same as the experimental possible energy ranges. 5,000 experiments 
with 500 total events on average in each experiment have been simulated. 

Remind that the measured energy spectra shown here are averaged over the simulated ex- 
periments. Five bins with linear increased bin widths have been used for binning generated 
signal and background events. As argued in Ref. [11], for reconstructing the one-dimensional 
WIMP velocity distribution function, this unusual, particular binning has been chosen in order 
to accumulate more events in high energy ranges and thus to reduce the statistical uncertainties 
in high velocity ranges. However, as shown in Sec. 2, for the determinations of ratios between 
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Figure 3: As in Figs. 2, except that an ^^''I target and the constant background spectrum have 
been used here. 

different WIMP couphngs/cross sections, one needs either events in the first energy bins or all 
events in the whole data sets. Hence, there is in practice no difference between using an equal 
bin width for all bins or a (linear) increased bin widths. 

It can be found in Figs. 1 that the shape of the WIMP scattering spectrum depends highly on 
the WIMP mass: for light WIMPs (m^ < 50 GeV), the recoil spectra drop sharply with increasing 
recoil energies, while for heavy WIMPs (m^ > 100 GeV), the spectra become flatter. In contrast, 
the exponential background spectra shown here depend only on the target mass and are rather 
flatter (sharper) for light (heavy) WIMP masses compared to the WIMP scattering spectra. This 
means that, once input WIMPs are light (heavy), background events would contribute relatively 
more to high (low) energy ranges, and, consequently, the measured energy spectra would mimic 
scattering spectra induced by heavier (lighter) WIMPs. Moreover, for heavy WIMP masses, 
since background events would contribute relatively more to low energy ranges, the estimated 
value of the measured recoil spectrum at the lowest experimental cut-off energy, r{Q^i^), could 
thus be (strongly) overestimated. 

Furthermore, in Figs. 2 we use a ^^F nucleus as detector material and fix the input WIMP 
mass as 100 GeV. Four different background ratios have been considered: no background (top 



left), 10% (top right), 20% (bottom left), and 40% (bottom right) background events in the 
analyzed data sets. It can be seen clearly that, for lighter nuclei e.g., F or Si, the WIMP 
scattering spectra are flatter than that of a Ge target, and more importantly, the exponential 
background spectra contribute (almost) only to low energy ranges (mostly into the first Q— bin 
and a bit into the second one) and does not affect the measured spectra (solid red histograms) 
in high energy ranges. In contrast. Figs. 3 show that, for heavier nuclei e.g., I or Xe, the WIMP 
scattering spectra are shaper than that of a Ge target and, since the constant background spectra 
contribute relatively mainly to high energy ranges (mostly into the last two Q— bin and a bit 
into that in the middle), the measured spectra in low energy ranges changes only very slightly. 
Gonsequently, the estimated value of r((5min) would only be (very) slightly overestimated. 

More detailed illustrations and discussions about the effects of residue background events 
with different spectrum forms on the measured energy spectrum and on the determination of 
the WIMP mass can be found in Ref. [17]. 

4 Results of the reconstructed ratios of WIMP— nucleon 
couplings/cross sections I: with exponential background 
spectra 

In this and the next sections I present simulation results of the reconstructed ratios between 
different WIMP-nucleon couplings/cross sections with mixed data sets from WIMP-induced 
and background events by means of the model-independent procedures described in Sec. 2.^^ 
Gonsidering the natural abundances of spin-sensitive detector materials (see Table 1), a ^^F and 
an ^^^I nuclei have been chosen as our targets. The threshold energies of all experiments have 
been assumed to be negligible^^ and the maximal cut-off energies are set the same as 100 keV. 
The exponential and constant background spectra given in Eqs. (38) and (37) have been used for 
generating background events in windows of the entire experimental possible ranges in this and 
the next sections, respectively. 2 (3) x 5,000 experiments have been simulated. Each experiment 
contains 50 total events on average. Note that all events recorded in our data sets are treated 
as WIMP signals in the analyses, although statistically we know that a fraction of these events 
could be backgrounds. 



4.1 Reconstructed (a^/ ap)^^ 

Gonsider at first the case of a dominant SD WIMP-nucleus interaction. Figs. 4 show the 
reconstructed fln/op ratios and the lower and upper bounds of their la statistical uncertainties 
estimated by Eqs. (13) and (A18) with n = —1 (dashed blue), 1 (solid red), and 2 (dash-dotted 
cyan) as functions of the input a^/a^ ratio. The background ratios shown here are no background 
(top left), 10% (top right), 20% (bottom left), and 40% (bottom right) background events in 
the analyzed data sets. Here the (minus)" solution has been used (as the "inner" solution) 
[10]. The mass of incident WIMPs has been set as 100 GeV. 

It can be found here that, firstly, the statistical uncertainty on the reconstructed a-^/a-p ratio 
with n = 2 is a little bit smaller than the uncertainty on the ratio reconstructed with 71 = 1, 

^^Note that, rather than the mean values, the (bounds on the) reconstructed ratios are always the median 
values of the simulated results. 

^^Different from our setup used in Refs. [9, 10], where the threshold energies have been set as 5 keV for all 
targets. 
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Figure 4: The reconstructed ttn/ttp ratios estimated by Eq. (13) and the lower and upper bounds 
of their la statistical uncertainties estimated by Eq. (A18) with n = — 1 (dashed blue), 1 (solid 
red), and 2 (dash-dotted cyan) as functions of the input On/ap ratio. Here the "— (minus)" 
solution has been used (as the "inner" solution) [10]. The mass of incident WIMPs has been set 
as 100 GeV and each experiment contains 50 total events on average. The other parameters are 
as in Figs. 2. and 3. 

and both of them are much smaller than that reconstructed with n = —1. Secondly, due to 
the non-negligible background ratio in the analyzed data sets, the reconstructed a^/a^ ratios 
become underestimated; the larger the background ratio the larger this systematic deviation of 
the reconstructed a^/a^. However, for the same data sets, the larger the n value (or, equivalently, 
the larger the moment of /i(f) used), the smaller this systematic deviation. This implies that the 
larger the background ratio, the more incompatibile between the an/op ratios reconstructed with 
different n. This (in) compatibility could thus offer us a simple check for the purity/availability 





Figure 5: As in Figs. 4, except that we estimate (an/ap)|^„ by Eq. (13) with the counting rates 
at the shifted points of the first Q— bin, r(^x,Y){Qs,i,{x,Y)) = r(x,Y),i- 

of our data sets. 

In Figs. 5 I show the reconstructed a^/a^ ratios and the lower and upper bounds of their la 
statistical uncertainties estimated by Eqs. (13) and (A18) with the counting rates at the shifted 
points of the first Q— bin, r(x,y)(Qs,i,(x,y)) = '"(X,y),i5 as functions of the input a^/a^ ratio^^. 
Different from the results shown in Figs. 4, the statistical uncertainty on the reconstructed a^/ Op 
ratio with n = — 1 is a little bit sma//er than those reconstructed with n = 2 and n = l.^'' But, 
the same as shown in Figs. 4, the larger the background ratio in the analyzed data sets, the 
more strongly the reconstructed On/ap ratios would be systematically underestimated, and the 

^^Labeled hereafter with an "sh" in the subscript. 

^^This is because we set simply the experimental threshold energies to be negligible in our simulations (see 
Ref. [10] for cases with non-negligible threshold energies). 





Figure 6: A comparison of the results shown in Figs. 4 estimated with r(x,y)(Qmin,(x,y) = 0) 
(dashed blue) and those shown in Figs. 5 with r(x,y)(Qs,i,(x,y)) (sohd red). Only the results 
estimated with n = 1 are shown here. 

larger the n value, the smaller this systematic deviation. However, the incompatibility between 
the fln/flp ratios reconstructed with different n is not so significant as shown in Figs. 4. 

As a comparison, I show the results shown in Figs. 4 estimated with r(^x,Y)iQmm,{x,Y) ~ 0) 
(dashed blue) and those shown in Figs. 5 with r(x,y)(Qs,i,(x,y)) (solid red) together (only cases 
with n = 1) in Figs. 6. It can be seen here clearly as well as from Figs. 4 and 5 that, because we 
set simply the experimental threshold energies to be negligible in our simulations, the statistical 
uncertainties on On/ap estimated with r(x,y)(Qs,i,(x,y)) are a little bit larger. However, as shown 
in Ref. [10], once threshold energies of analyzed data sets are non-negligible, the statistical 
uncertainties on Cn/ap estimated with r(x,y)(Qs,i,(x,y)) could be ~ 16% smaller than those 
estimated with r(^x,Y){Qmin,{x,Y) = 5 keV). Moreover, as shown here as well as in Figs. 4 and 
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Figure 7: The reconstructed Cn/ap ratios estimated by Eq. (13) and the lower and upper bounds 
of their la statistical uncertainties estimated by Eq. (A18) with n = —1 (dashed blue), 1 (solid 
red), and 2 (dash-dotted cyan) as functions of the input WIMP mass m^. The input a^/ap ratio 
has been set as 0.7. The other parameters are as in Figs. 4. 

5, for the cases with r(^x,Y)iQs,i,{x,Y)) systematic deviations caused by residue background 
events are much smaller. 

Quantitatively, Figs. 4 to 6 show that, with even ~ 20% - 40% residue background events in 
the analyzed data sets, one could in principle still reconstruct the ratio between the SD WIMP 
couplings on neutrons and that on protons pretty well; for a WIMP mass = 100 GeV and 
dn/dp = 0.75, by using Eq. (13) with r(x,y)(Qmin,(x,y) = 0) (with r^x,Y)iQs,i,{x,Y))) and = 1 to 
analyze data sets of a 10% (20%) background ratio, the systematic deviation could be < — 12% 
( < - 14%) with an ~ 55% (~ 60%) statistical uncertainty^^. 

^^In Refs. [9, 10] another combination of detector materials: ^^Ge + "^^Cl has also been considered. With 
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Figure 8: As in Figs. 7, except that we estimate (an/ap)|^„ by Eq. (13) with the counting rates 
at the shifted points of the first Q— bin, r(^x,Y){Qs,i,{x,Y)) = r(x,Y),i- 



The expression (13) for estimating the ratio between two SD WIMP-nucleon couphngs is 
independent of the WIMP mass. However, as discussed in Sec. 3 and in Ref. [17] in more detail, 
an exponential-hke residue background spectrum could cause an over-/underestimate of the 
reconstructed WIMP mass for light/heavy WIMPs. Hence, in order to check the WIMP-mass 
independence of the reconstructed ctn/op ratio with non-negligible background events, in Figs. 7 
and 8 I show the reconstructed results with n = —1 (dashed blue), 1 (solid red), and 2 (dash- 
dotted cyan) as functions of the input WIMP mass m^. The input a^/ap ratio has been set as 
0.7. 



the same setup used here, expect non-negUgible experimental threshold energies (Qmin — 5 keV), the statistical 
uncertainties on the a-n/a-p ratios reconstructed with Gc + CI are only ^ 1/5 — 1/3 of those reconstructed with 
F + I. For detailed discussions about the reasons of this difference, see Ref. [10]. 
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Figure 9: A comparison of the results shown in Figs. 7 estimated with r(x,y)(Qmin,(x,y) = 0) 
(dashed blue) and those shown in Figs. 8 with r(^x,Y)iQs,i,(x,Y)) (sohd red). Only the results 
estimated with n = 1 are shown here. 

It can be seen here that, for WIMP masses > 50 GeV, except the statistical uncer- 
tainty estimated with r(x,y)(Qmm,(x,y)) and n = —1 (dashed blue curves labeled as (an/ap)r°, _i 
in Figs. 7), the reconstructed a^/ap ratios as well as their statistical uncertainties are indeed 
(almost) independent of the WIMP mass. However, as discussed above, the larger the back- 
ground ratios in our data sets, the more strongly underestimated the reconstructed a^/ap ratios 
for all input WIMP masses. And, as shown in Figs. 9, the a^/ap ratios reconstructed with 
f'{X,Y) (Qmin, ,{x,Y)) (dashed blue) would be more strongly underestimated as those reconstructed 
with r(x,y)(Qs,i,(x,y)) (solid red). Nevertheless, with data sets of < 20% residue background 
events, the reconstructed la statistical uncertainty intervals could in principle always cover the 
input (true) On/ap ratios pretty well. 





Figure 10: As in Figs. 6, except that the SI WIMP-nucleon cross section and the SD WIMP- 
proton couphng have been set as cr^p = 10~^ pb and Op = 0.1, respectively. 

On the other hand, for WIMP masses < 50 GeV, interestingly. Figs. 7 show that the On/ ctp 
ratio reconstructed with r(x,Y)(Qmin,(x,y)) and n = —1 becomes the best result with an < 20% 
background ratio the la statistical uncertainty interval could still cover the input (true) a^/a-^ 
ratios well. In contrast, all other five reconstructed an/op ratios shown in Figs. 7 and 8 are 
(strongly) underestimated. 




Figure 11: The reconstructed an/cip ratios estimated by Eqs. (13) (dashed blue, n = 1) and (28) 
(sohd red) and the lower and upper bounds of their la statistical uncertainties estimated by 
Eqs. (A18) and (A22) with r(^x,Y,z){Qs,i,(x,Y,z)) as functions of the input a^/ap ratio. Besides 
^^F and ^^'^I, ^^Si has been chosen as the third target for estimating Cp^(x,y) by Eqs. (29a) and 
(29b). The other parameters are as in Figs. 10. 

4.2 Reconstructed (an/ap)^^+^^ 

In this and the next subsections I consider the case with a non-negligible SI WIMP-nucleus 
interaction. The input SI WIMP-nucleon cross section and the input SD WIMP-proton coupling 
have been set as = 10~^ pb and Op = 0.1, respectively. 

"'^^Remind that this conclusion holds only when the experimental threshold energies of the analyzed data sets 
can be negligible. Otherwise, as shown in Ref. [10], the an/op ratio reconstructed with f(x.v') (Qmin.(x,F)) 
71 = —1 would also be strongly underestimated. 




m^jn [GeV] 



500 700 1000 



70 1 00 200 



500 700 1000 



I, Q^g, < 100 keV, 0^3^ bg """^ ^^'^< 2 X 50 events (20% exponential bg), c^p= 10"" pb. = 0.1, a^/Sp = 0.7 



<:100 keV, 2x50 





1 [GeV] 



500 700 1000 



70 1 00 200 

m i„ [GeV] 



500 700 1000 



Figure 12: As in Figs. 9, except that the SI WIMP-nucleon cross section and the SD WIMP- 
proton couphng have been set as cr^p = 10~^ pb and Op = 0.1, respectively. 

At first I show in Figs. 10 the reconstructed a^/a^ ratios estimated by Eq. (13) and the 
lower and upper bounds of their la statistical uncertainties estimated by Eq. (A18) with 
'^(x,r)(Qmin.(x,y)) (dashed blue) and with r(x,y)(Qs,i,(x,r)) (sohd red) as functions of the in- 
put On/op ratio together (only cases with n = 1). Since for this simulation setup the SD 
WIMP-nucleus cross section doesn't really dominate over the SI one, by using data sets with 
pure WIMP signals (no background events, top left frame), the reconstructed On/ap ratios are 
(strongly) overestimated, especially for the input On/op < 0. However, once some background 
events mix into our data sets, the reconstructed a^/a^ ratios become smaller (and even closer to 
the input (true) values, cf. Figs. 4 to 6). This seems to indicate that, with data sets of < 40% 
residue background events in the analyzed data sets, one could in principle reconstruct the ratio 
between two SD WIMP-nucleon couplings by simply assuming a dominant SD WIMP-nucleus 




Figure 13: The reconstructed a-a/a^ ratios estimated by Eq. (28) and the lower and upper bounds 
of their la statistical uncertainties estimated by Eq. (A22) with r(x,y,z)(Qmm,(x,y,z)) (dashed 
blue) and with r(x,y,z)(Qs,i,(x,y,z)) (sohd red) as functions of the input WIMP mass m^. Besides 
^^F and ^^''^I, ^^Si has been chosen as the third target for estimating c^^{x,y) by Eqs. (29a) and 
(29b). The other parameters are as in Figs. 12. 

interaction (even though this assumption is not correct) and using Eq. (13). 

Nevertheless, once we consider both SI and SD WIMP interactions and thus use Eqs. (28) and 
(A22) to analyze the same data sets, as shown in Figs. 11, the an/op ratios could be reconstructed 
(much) better with data sets of < 20% residue background events, even though the statisti- 
cal uncertainties are pretty large for input a^Ja^ < 0. For a WIMP mass = 100 GeV and 
On/«p = 0.75, by using Eq. (28) with r(x,y,z)(<5s,i,(x,y,z)) to analyze data sets with a 10% back- 
ground ratio, the systematic deviation could be ~ —12% with an ~ 50% statistical uncertainty. 
Moreover, with an increased background ratio, the (in) compatibility between the reconstructed 




Figure 14: As in Figs. 11, except that the SI WIMP-nucleon cross section has been set as 



all = 10-^° pb. 



ttn/ttp ratios estimated by Eqs. (13) (dashed blue, n = 1) and (28) (sohd red) becomes larger. 
Hence, as mentioned in the previous subsection, one could/should compare results reconstructed 
under different assumptions, with both r(^x,Y,z)iQmm,(x,Y,z)) and r(_^x,Y,z){Qs,i,{x,Y,z)) with 
different n values, in order to check the purity of the analyzed data sets (as well as the dominance 
of the SI or the SD WIMP interaction [10]). 

On the other hand, as done in the previous subsection, in order to check the WIMP-mass in- 
dependence of the reconstructed results with non-negligible background events, I show the a^/ Op 
ratios estimated by Eq. (13) with r(x,y)(Qmm,(x,y)) (dashed blue) and with r(x,y)(Qs,i,(x,y)) 
(solid red) as functions of the input WIMP mass together (only cases with n = 1) in 
Figs. 12 as well as those estimated by Eq. (28) with r(x,y,z)(Qmm,(x,y,z)) (dashed blue) and with 
''^{x,Y,z){Qs,i,{x,Y,z)) (solid red) as functions of the input WIMP mass in Figs. 13. 





Figure 15: The reconstructed a^/a^^ ratios and the lower and upper bounds of their la sta- 
tistical uncertainties as functions of the input On/flp ratio. The dashed blue curves indicate the 
ratios estimated by Eq. (26) with On/ap estimated by Eq. (28) {not by Eq. (13)), whereas the 
solid red curves indicate the ratios estimated by Eq. (31). ™Ge has been chosen as the second 
target having only the SI interaction with WIMPs and combined with ^'^Na for using Eq. (31). 
Parameters are as in Figs. 11. Note that the input a^/ap ratio ranges only between —1 and 2. 

As shown earlier, for WIMP masses ^ 50 GeV, the reconstructed a^/ Op ratios are indeed 
(almost) independent of m^. And, with data sets of < 20% background ratios, using Eq. (28) 
(solid red in Figs. 13) and estimating with r(^x,Y,z){Qs,i,{x,Y,z)) could offer the best reconstructed 
On/ap ratios. Meanwhile, for WIMP masses < 50 GeV, Figs. 13 show also that the an/op 
ratios reconstructed by Eq. (28) with r(x,y,z)(Qmin,(x,y,z)) could be the best results^° (with 
however a pretty large statistical uncertainty): for a WIMP mass = 25 GeV and an input 



See footnote 19. 
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Figure 16: The reconstructed a^/a^^ ratios and the lower and upper bounds of their la sta- 
tistical uncertainties as functions of the input WIMP mass m^. The input an/a^ = 0.7, the 
other parameters and notations are as in Figs. 15. Note that the input WIMP mass starts from 
25 GeV. 

fln/ctp = 0.7, by using Eq. (28) with r(x,y,z)(Qmm,(x,y,z)) to analyze data sets of a 10% background 
ratio, the systematic deviation could be ~ —25% with an ~ 88% statistical uncertainty. 

Furthermore, in Figs. 14 we reduce the SI WIMP-nucleon cross section two orders of magni- 
tude smaller: cx^p = 10~^° pb. It can be found here that, although the a^/a^ ratios reconstructed 
by Eqs. (13) (dashed blue, n = 1) and (28) (sohd red) with r(x,y,z)(Qs,i,(x,y,z)) niatch much 
better than those shown in Figs. 11, with an increased background ratio the a^/ap ratios could 
be a bit more strongly underestimated by using Eq. (28). Nevertheless, Figs. 11, Figs. 13, and 
Figs. 14 show that, it doesn't matter whether the SD WIMP-nucleus interaction really domi- 
nates over the SI one or not, by using Eq. (28) one could in principle always reconstruct the 




Figure 17: The reconstructed a^/a^^ ratios and the lower and upper bounds of their la sta- 
tistical uncertainties as functions of the input a^/ap ratio. The dashed blue curves indicate the 
ratios estimated by Eq. (26) with On/ctp estimated by Eq. (28), whereas the solid red curves 
indicate the ratios estimated by Eq. (31). ^^Ge has been chosen as the second target having only 
the SI interaction with WIMPs and combined with ^'^^Xe for using Eq. (31). Parameters are as 
in Figs. 15. Note that the input an/op ratio ranges only between —1 and 2. 

ratio of the SD WIMP coupling on neutrons to that on protons with data sets of < 20% residue 
background events pretty well. But, the larger the relative strength between the SD WIMP- 
nucleus interaction to the SI one, the smaller the systematic deviations as well as the statistical 
uncertainties. For = 10~^° pb with a WIMP mass = 100 GeV and a^/a^ = 0.75, by 
using data sets of a 10% background ratio, the systematic deviation could be ~ —10% with an 
~ 43% statistical uncertainty. 




Q^ax < keV, Qmax, bg < ^^^V, 3/2 X 50 events (20% exponential bg), o^p= 10"" pb, Bp = 0.1, a^/^p = 0-'' 



AMIDAS http://pisrvO.pit.physik.utii-tuebingen.de/darkmattet/amidas/ 



in / rv 1 



50 70 100 200 300 

[GeV] 



500 700 1 000 



Qn,aj, < 100 keV, Q^a^ bg 100 l^^^, 3/2 X 50 events (40% exponential bg), o^p= 10'" pb, ap = 0.1, a^/ap = 0.7 



. ANIIDAS htlp;//pistv0.pil.physik. uni-tuBbingen.de/dathmatter/amidas/ 



50 70 100 200 300 

m i„ [GeV] 



500 700 1000 



Figure 18: The reconstructed a^/a^^ ratios and the lower and upper bounds of their la sta- 
tistical uncertainties as functions of the input WIMP mass m^. The input an/a^ = 0.7, the 
other parameters and notations are as in Figs. 17. Note that the input WIMP mass starts from 
25 GeV. 

4.3 Reconstructed cr^(Jj,n)/'^xp 

In Figs. 15 and 16 I show the reconstructed a^/a^^ ratios and the lower and upper bounds 
of their la statistical uncertainties as functions of the input a^/a^ ratio as well as of the input 
WIMP mass m^, respectively. The dashed blue curves indicate the ratios estimated by Eq. (26) 
with an/ctp estimated by Eq. (28) {not by Eq. (13)), whereas the solid red curves indicate the 
ratios estimated by Eq. (31). ^^Ge has been chosen as the second target having only the SI 
interaction with WIMPs and combined with ^^Na for using Eq. (31). 

It can be seen here that, interestingly, while the a^/a^^ ratios reconstructed with F + 1 + Si 
targets become as usual more and more strongly underestimated with an increased background 




Figure 19: As in Figs. 5, except that the constant background spectrum has been used here. 

ratio, those reconstructed with Na + Ge targets become in contrast more and more strongly over- 
estimated. Nevertheless, for WIMP masses > 25 GeV, with < 20% residue background 
events, the overlap of the la statistical uncertainty intervals estimated by two target combina- 
tions could cover the input (true) cr^/c^n I'atio well. For a WIMP mass = 100 GeV and 
«n/«p = 0.7 (the theoretical value of o-|p/c^fp = 3.49 x 10^), by using F + I + Si (Na + Ge) targets 
with data sets of a 20% background ratio, one could in principle reconstruct the a^/a^ ratio 
with an ~ —20% (~ +21%) systematic deviation and an ~ 31% (~ 32%) statistical uncertainty. 

On the other hand, in Figs. 17 and 18 I show the reconstructed a^/a^^ ratios and the lower 
and upper bounds of their la statistical uncertainties as functions of the input an/op ratio as 
well as of the input WIMP mass m^, respectively. The dashed blue curves indicate the ratios 
estimated by Eq. (26) with a^/ap estimated by Eq. (28), whereas the sohd red curves indicate 
the ratios estimated by Eq. (31). ''^Ge has been chosen as the second target having only the SI 
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Figure 20: As in Figs. 7, except that the constant background spectrum has been used here. 



interaction with WIMPs and combined with ^^^Xe for using Eq. (31). 

It can be found here that, more interestingly, while the ratios reconstructed with 

F -f I + Si targets become more and more strongly underestimated with an increased background 
ratio for all input a^/ap values, those reconstructed with Xe + Ge targets become more and more 
strongly underestimated for a^/ap > and more and more strongly overestimated for a^/a^ < 0. 
Nevertheless, for WIMP masses > 25 GeV, with < 20% residue background events, the 
(overlap of the) la statistical uncertainty intervals estimated by two target combinations could 
cover the input (true) a^^/a^l^ ratio well. For a WIMP mass 



100 GeV and a^/ap = 0.7 
(the theoretical value of Cyn/'^vp — ^■'^^ ^ 10^), by using F + I + Si (Xe + Ge) targets with data 



sets of a 10% background ratio, one could in principle reconstruct the ratio with an 

~ —32% (~ —3.2%) systematic deviation and an ~ 120% (~ 62%) statistical uncertainty. 




Figure 21: As in Figs. 8, except that the constant background spectrum has been used here. 

5 Results of the reconstructed ratios of WIMP— nucleon 
couplings/cross sections II: with constant background 
spectra 

In this section I show simulation results with residue background events generated by the con- 
stant spectrum given in Eq. (37) and compare them with those shown in the previous section. 
Some general rules about the effects of different background sources will also be discussed. 




Figure 22: As in Figs. 11, except that the constant background spectrum has been used here. 
5.1 Reconstructed (a^/ cip)^^ 

As in the previous section, I consider at first the case of a dominant SD WIMP-nucleus inter- 
action. 

In Figs. 19 I show the reconstructed On/op ratios and the lower and upper bounds of their la 
statistical uncertainties estimated by Eqs. (13) and (A18) with n = —1 (dashed blue), 1 (solid 
red), and 2 (dash-dotted cyan) as functions of the input a^/ttp ratio. Note that only results with 
the counting rates at the shifted points of the first Q— bin, r(x,y)(Qs,i,(x,y)), are shown here, 
because their systematic deviations due to the non-negligible background events as well as their 
statistical uncertainties (for non-zero experimental threshold energies) are (much) smaller^^. 

It can be found here that, firstly, in contrast to the results shown in Figs. 4 to 6, the a^/ap 



'^^See Figs. 6 and discussions there. 





Figure 23: As in Figs. 12 and 13, except that the constant background spectrum has been used 
here. The dashed blue (sohd red) curves indicate the reconstructed an/ap ratios estimated by 
Eq. (13) with n = 1 (Eq. (28)) with r(x,Y)(<5mm,(x,y)) as functions of the input On/flp ratio. 

ratios reconstructed with n = 1 (sohd red) and 2 (dash-dotted cyan) are now overestimated. 
This should be caused by the background contribution to high energy ranges. Remind that, as 
shown in Figs. 2 and 3, while an exponential(-like) background spectrum contributes (almost) 
only/mainly to low energy ranges, a constant/(approximately) flat one contributes mainly to 
high energy ranges. This indicates in turn that, while the reconstructed a^/a^ ratio would be 
underestimated by using data sets with residue background events existing mostly in low energy 
ranges, e.g., electronic noise or incompletely charged surface events, the On/ap ratio would be 
overestimated hj using data sets with backgrounds relatively mainly in high energy ranges, e.g., 
cosmic rays and cosmic-ray induced 7-rays with energies of (9(100) keV. 

However, interestingly and importantly, the a^/a^ ratio reconstructed with n = —1 (dashed 
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Figure 24: As in Figs. 15, except that the constant background spectrum has been used here. 

blue) in Figs. 19 is almost not affected by the constant background spectrum! The reason can 
be understood as follows. The (an/«p)^°_i sh ratio has been estimated with 'R-j-i^{xx) given in 
Eq. (21) and is in fact a function of only r (^x,y){Q s,i,{x,y)) ■ Since the constant background spectra 
contribute only (very) small amounts to the lowest energy ranges (see Figs. 3), r(^x,Y){Qs,i,{x.Y)) 
estimated by using events recorded in the first Q— bins would not change or only very slightly, 
and in turn also the reconstructed a^/a^ ratio. 

In Figs. 20 and 21, I show the reconstructed a^/a^ ratios as functions of the input WIMP 
mass. Both of them show that, while the a^/a^ ratios reconstructed with n = 1 (solid red) and 
2 (dash-dotted cyan) are more and more strongly overestimated with an increased background 
ratio, the ratios reconstructed with n = —1 (dashed blue) just become a little bit smaller (~ 15% 
smaller for a background ratio of 40% and an input WIMP mass of 1 TeV, compared to the results 
with no background events) and the statistical uncertainties also grow only slightly ( ^ 45% 





Figure 25: As in Figs. 16, except that the constant background spectrum has been used here. 



larger) . 

5.2 Reconstructed (an/ap)^^+^° 

In this subsection I consider the reconstruction of the ctn/op ratio with a non-zero SI WIMP- 
nucleus cross section. 

As shown in the previous subsection, Figs. 22 and 23 show that, while with an increased 
background ratio the an/flp ratios reconstructed by Eq. (13) (dashed blue, with n = 1 and 
'^(x,y)(<5min,(x,y) = 0) or r(x,y)(Qs,i,(x,y))) are more and more strongly overestimated caused by 
contributions of the constant background spectrum to high energy ranges, those reconstructed 
by Eq. (28) (solid red) just become a little bit smaller (~ 9% smaller for a background ratio of 
40% and an input WIMP mass of 100 GeV, compared to the results with no background events) 





Figure 26: As in Figs. 17, except that the constant background spectrum has been used here. 

and the statistical uncertainties also grow only slightly (~ 34% larger). 

This is simply because that {a^/ ap)^~^^^ given in Eq. (28) depends only on Cp_(x,y), which are 
in turn just the functions of ?"(x,y,z)(Qmin,(x,y,z)) or r(^x,Y,z){Qs,i,{x,Y,z)) and can be estimated by 
using events in the lowest energy ranges. Hence, once residue background events exist regularly 
between the experimental minimal and maximal cut-off energies or (even better) (mostly) in 
high energy ranges, one can in principle estimate r(x,y,z)(Qmm,(x,y,z)) and r(x,y,z)(Qs,i,(x,y,z)) as 
well as the ratio between two SD WIMP-nucleon couplings (pretty) precisely by using Eq. (28) 
without worrying about the non-negligible backgrounds. 





Figure 27: As in Figs. 18, except that the constant background spectrum has been used here. 
5.3 Reconstructed <t^P x/ct^^ 

Finally, we check the reconstruction of the ratios between the SD and SI WIMP-proton (neutron) 
cross sections by using Eqs. (26) and (31) with the constant background spectrum. 

In Figs. 24 and 25 we can find that, as observed in Sec. 4.3, while the ratios re- 

constructed with F + I + Si targets by Eq. (26) with fln/^p reconstructed by Eq. (28) (dashed 
blue) become more and more strongly underestimated with an increased background ratio, those 
reconstructed with Na + Ge targets by Eq. (31) (solid red) become more and more strongly over- 
estimated. Meanwhile, Figs. 26 and 27 show that, while the a^/a^^ ratios reconstructed with 
F + I + Si targets (dashed blue) become more and more strongly underestimated with an in- 
creased background ratio for all a-^/a-p values, those reconstructed with Xe + Ge targets (solid 
red) become more and more strongly underestimated for a^/ap > and more and more strongly 



overestimated for an/op < 0. 

However, the shifts caused by events with the constant background spectrum are much 
smaller than those caused with the exponential one. Quantitatively, for an input WIMP mass 
of 100 GeV, an input SI WIMP-nucleon cross section of 10^*^ pb, an input SD WIMP-proton 
coupling ttp = 0.1, and an input a-^/ap = 0.7, by using 2 (3) data sets of a 40% background 
ratio, the a^/a^^ ratios would be reconstructed 4% larger (7% smaller) with < 20% larger 
statistical uncertainties, whereas the o"^n/c^p ratios would be reconstructed 14% (30%) smaller 
with ~ 30% larger statistical uncertainties, compared to the results with no background events. 

These results indicate that, once residue background events exist regularly between the ex- 
perimental minimal and maximal cut-off energies or (mostly) in high energy ranges and one can 
therefore in principle estimate r^x,Y,z){Qniin,(x,Y,z)) and r(^x,Y,z){Qs,i,{x,Y,z)) (pretty) well, the 
ratios between the SD and SI WIMP-nucleon cross sections could then be estimated (pretty) 
precisely without worrying about the non-negligible backgrounds. 

6 Summary and conclusions 

In this paper I reexamine the model-independent data analysis methods introduced in Rcfs. [9, 
10] for the determinations of ratios between different WIMP-nucleon couphngs/cross sections 
from data (measured recoil energies) of direct Dark Matter detection experiments directly by 
taking into account a fraction of residue background events, which pass all discrimination criteria 
and then mix with other real WIMP-induced events in the analyzed data sets. These methods 
require neither prior knowledge about the WIMP scattering and different possible background 
spectra nor about the WIMP mass; the unique needed information is the recoil energies recorded 
in two (or more) direct detection experiments. 

I considered at first the case of a dominant SD WIMP-nucleus interaction. Our simula- 
tions show that, due to the contribution of non-negligible residue background events in the 
analyzed data sets to low/high energy ranges, the reconstructed a^/ap ratios would be under- 
/overestimated; the larger the background ratio the larger these systematic deviations of the 
reconstructed On/op. But, by estimating the counting rates at the shifted points, instead of at 
the experimental minimal cut-off energies, one could (strongly) alleviate these systematic devi- 
ations as well as reduce the statistical uncertainties (with non-negligible experimental threshold 
energies). By using data sets of ~ 20% - 40% residue background events, the ratio between 
the SD WIMP couphngs on neutrons and on protons could in principle still be reconstructed 
pretty well: for a WIMP mass — 100 GeV and On/op = 0.75, with data sets of a 10% (20%) 
background ratio, the systematic deviation could be < — 12% ( < — 14%) with an ~ 55% 
(~ 60%) statistical uncertainty. 

Then I turned to consider the general combination of the SD WIMP-nucleus interaction with 
a non-negligible SI one. Our simulations show that, by combining three (two spin-sensitive) 
target nuclei, our method can be used to reconstruct the ratio of the SD WIMP couphng on 
neutrons to that on protons with data sets of < 20% residue background events. And, more 
importantly, it doesn't matter whether the SD WIMP-nucleus cross section really dominates over 
the SI one or not. But, the larger the relative strength between the SD WIMP-nucleus interaction 
to the SI one, the smaller the systematic deviations as well as the statistical uncertainties. 

I considered also the reconstruction of the ratios between the SD and SI WIMP-nucleon cross 
sections. For WIMP masses > 25 GeV, by using either the reconstructed a^^/ap ratio or a 
combination of a spin-sensitive nucleus with an only Sl-sensitive one with data sets of < 20% 
background ratios, while the cr^/o'^ ratio could be reconstructed with a ~ 30% statistical 



uncertainty, one could (only) estimate the order of magnitude of the / cr^p ratio (because of 
an ~ 120% statistical uncertainty). 

Moreover, our simulations show also the WIMP-mass independence of the reconstructed 
dn/dp ratios with non-negligible background events, especially for WIMP masses > 50 GeV. 
For WIMP masses < 50 GeV, the Cn/op ratios could be (strongly) underestimated, even with zero 
background events. However, this underestimate could be alleviated/corrected by decrease the 
experimental minimal cut-off energies of the analyzed date sets (to be negligible). Then, with 
data sets of < 20% residue background events, one could still reconstruct the fln/o-p ratios pretty 
well, by either assuming a dominant SD WIMP interaction or using the general combination of 
the SI and SD cross sections. But, the statistical uncertainty could be pretty large, once the SD 
WIMP-nucleus interaction doesn't dominate over the SI one. 

Furthermore, it has also been found that, firstly, by taking different assumptions about the 
relative strength between the SI and SD WIMP-nucleus interactions, and/or using different 
moments of the one-dimensional WIMP velocity distribution function, at either the experimen- 
tal minimal cut-off energies or the shifted energy points, there would be an (in) compatibility 
between different reconstructed On/op ratios; with an increased background ratio, the incom- 
patibility between the reconstructed results would become larger. Hence, this (in) compatibility 
could allow us to check the purity/availability of the analyzed data sets (as well as the dominance 
of the SI or SD WIMP interaction). 

Secondly and more importantly, our simulations with a constant background spectrum indi- 
cate that, once residue background events exist regularly between the experimental minimal and 
maximal cut-off energies or (even better) (mostly) in high energy ranges, one could in principle 
estimate the counting rates of the recoil spectrum of only WIMP-induccd events at the exper- 
imental threshold energies (pretty) precisely. Then the ratio between two SD WIMP-nucleon 
couphngs as well as the ratios between the SD and SI WIMP-nucleon cross sections could be 
estimated (pretty) precisely without worrying about the non-negligible backgrounds. 

In summary, as the forth part of the study of the effects of residue background events in 
direct Dark Matter detection experiments, we considered the determinations of ratios between 
different WIMP-nucleon couplings/cross sections. Our results show that, with currently running 
and projected experiments using detectors with 10~^ to lO"-*^^ pb sensitivities [32, 14, 33, 34] 
and < 10~^ background rejection ability [13, 15, 16, 12], once two or more experiments with 
different spin-sensitive target nuclei could accumulate a few tens events (in one experiment), 
we could in principle already estimate the relative strengths of couplings/cross sections of Dark 
Matter particles on ordinary matter with a reasonable precession, even though there could be 
some background events mixed in our data sets for analyses. Moreover, although two forms for 
background spectrum considered in this work is rather naive, the nuclear form factors for the SD 
WIMP interaction with different target nuclei are also more complicated as the simple thin-shell 
form used in our simulations, and the relative signs of the (ratios of the) expected/measured 
proton/neutron group spins of the used target nuclei could also change the reconstructed results 
(to be larger or smaller, underestimated or overestimated), one should be able to extend our 
observations/discussions to predict the effects of possible background events in their own ex- 
periments. Hopefully, this will not only encourage our experimental colleagues to present their 
(future) results in the parameter space of Dark Matter particles, but also help them to check 
the purity of their data sets, to understand (residue) background events in their experiments, as 
well as to improve their background discrimination techniques. 
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A Formulae needed in Sec. 2 

Here 1 list all formulae needed for the model-independent data analysis procedures used in Sec. 2. 
Detailed derivations and discussions can be found in Refs. [11, 10]. 



A.l Estimating r(Qmin) and /n(Qmin, Qmax) 

First, consider experimental data described by 



O — hrt <_ O < O +^ 



1, 2, 



, Ar„, n= 1, 2, 



B. 



(Al) 



Here the total energy range between Qmin and Qmax has been divided into B bins with central 
points Qn and widths 6„. In each bin, Nn events will be recorded. Since the recoil spectrum 

dR/dQ is expected to be approximately exponential, the following ansatz for the measured ieco][ 
spectrum [before normalized by the experimental exposure £) in the nth bin has been introduced 
[11]: 
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Here r„ is the standard estimator for (dR/ dQ)e^pt Q — Qn- 
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kn is the logarithmic slope of the recoil spectrum in the nth Q— bin, which can be computed 
numerically from the average value of the measured recoil energies in this bin: 
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The error on the logarithmic slope kn can be estimated from Eq. (A4) directly as 
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Qs,n in the ansatz (A2) is the shifted point at which the leading systematic error due to the 
ansatz is minimal [11], 



Qs,n Qn I, 



(A8) 



Note that Qs,n differs from the central point of the nth bin, Q^. Prom the ansatz (A2), the 
counting rate at Q — Qmm can be calculated by 



and its statistical error can be expressed as 
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Finally, since all are determined from the same data, they are correlated with 

iVtot Q{n+m~2)/2 

cov(4,/J = E^^, 



(All) 



(A12) 



where the sum runs over all events with recoil energy between Qmin and Qmax- And the correlation 
between the errors on r{Q-[am), which is calculated entirely from the events in the first bin, and 
on /„ is given by 
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note that the sums here only count in the first bin, which ends at Q = Qmin + ^i- 

On the other hand, with a functional form of the recoil spectrum (e.g., fitted to experimental 
data), {dR/ dQ)e:x.pt, one can use the following integral forms to replace the summations given 
above. Firstly, the average Q— value in the nth bin defined in Eq. (A5) can be calculated by 
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For /n(Qmm, Qmax) givcn in Eq. (11), we have 

/•Qmax g("-i)/2 /dR\ 

-^n(Qmin) Qmax) — I ,-,9 dQ , 

and similarly for the covariance matrix for 7„ in Eq. (A12), 



(A14) 



rQm. 

COv(/„, Im) = / 

m in 



g("+'"-^)/^ /dR\ 

F\Q) [dQj 



dQ. 



(A15) 



(A16) 



expt 



Remind that (di?/d(5)expt is the measured recoil spectrum before normahzed by the exposure. 
Finally, /^(Qmin, Qmin + ^'i) needed in Eq. (A13) can be calculated by 
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Note that, firstly, r((5mm) and /„((5min, <5min + ^'i) should be estimated by Eqs. (A9) and (A17) 
with ri, ki and Qs,i estimated by Eqs. (A3), (A4), and (A8) in order to use the other formulae for 
estimating the (correlations between the) statistical errors without any modification. Secondly, 
r((5min) and /n((5min, Qmax) estimated from a scattering spectrum fitted to experimental data 
are usually not model-independent any more. Moreover, for the use of Eqs. (11), (A12), (A15), 
(A16), and (A17) the elastic nuclear form factor F'^{Q) should be understood to be chosen for 
the SI and SD WIMP-nucleon cross section correspondingly. 



A. 2 Statistical uncertainty on {a^/a^)^ 
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By using the standard Gaussian error propagation, the statistical uncertainty on {an/ap)j^^ 
estimated by Eq. (13) can be expressed as 
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Here a short-hand notation for the six quantities on which the estimate of {cin/(ip)±^n depends 
has been introduced: 

Cl,X = In,X , C2,X = h,X , C3,x = rx(Qmin,x) ] (Al9) 

and similarly for the Q^y. Estimators for cov(Q,Cj) have been given in Eqs. (A12) and (A13). 
Explicit expressions for the derivatives of 7ln,x given in Eq. (17) with respect to Ci^x are: 
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(A20a) 
(A20b) 



7^„,x ; (A20c) 



explicit expressions for the derivatives of TZn,Y with respect to Cj,y can be given analogously. Note 

that, firstly, factors TZn,{x,Y) appear in all these expressions, which can practically be cancelled 
by the prefactors in the bracket in Eq. (A18). Secondly, all the /o,{x,y) and In,(x,Y) should be 
understood to be computed according to Eq. (11) or (A15) with integration limits Qmm and 
Qmax specific for that target. 

Similarly, expressions for the derivatives of TZc^^x can be computed from Eq. (16) as 
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(A21b) 



and similarly for the derivatives of 7la,Y- Remind that factors appearing here can also 

be cancelled by the prefactors in the bracket in Eq. (A18). 

A. 3 Statistical uncertainty on (an/«p)±''"^^ 

Prom the expression (28), the statistical uncertainty on (an/ap)|^+^° can be given by 
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Here, from the first and second lines of the expression (28), we have, 
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(for Sn/p,X < Sn/p,Y), 
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(for Sn/p,x < Sn/px)- 



Then, from the definition (29a) of Cp^x, one can get directly 
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Similarly, from the definition (29b) of Cp_y, we have 
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A. 4 Statistical uncertainty on cr^/cr^ 



Since Cp and defined in Eqs. (23) and (32) are functions of Un/ap, once the an/op ratio has 
been estimated (from e.g., some other direct detection experiments by Eq. (13) by assuming a 
dominant SD WIMP-nucleus interaction), c'^p/c^p can then be estimated by Eq. (26) with the 
following statistical uncertainty: 
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where 
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Here, from the expression (26) for estimating a^^/a^l^, its derivatives with respect 
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to Cp,(x,y) 

(A28a) 
(A28b) 



And the derivatives of cr^p/c^p ^^^^ respect to 7?.m,(x,y) are 
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Meanwhile, from expression (23) for Cp one can find that 
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and, since we estimate in fact always a^/ttp, one needs practically 
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On the other hand, for estimating a-^/ap by Eq. (28), the statistical uncertainty on 
can be expressed as 
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with 9((7^°/(7^^)/9(an/ap) given in Eq. (A27) and 
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for i — X, Y, Z. Moreover, by using Eq. (33) to eliminate the On/op dependence of a^/a^^, 
the statistical uncertainty given in Eq. (A32) can be reduced to 
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